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ABSTRACT 

We present the first large-scale study of E+A galaxies that incorporates photometry 
in the ultraviolet {UV) wavelengths. E+A galaxies are 'post-starburst' systems, with 
i-S^ . strong Balmer absorption lines indicating significant recent star formation, but without 

Dh' [Oil] and Ha emission lines which are characteristic of ongoing star formation. The 

starburst that creates the E+A galaxy typically takes place within the last Gyr and 
^ • creates a high fraction (20-60 percent) of the stellar mass in the remnant over a short 

^ ' timescale (< 0.1 Gyrs). We find a tight correlation between the luminosity of our 

, E+A galaxies and the implied star formation rate (SFR) during the starburst. While 

low-luminosity E+As (M(z) > -20) exhibit implied SFRs of less than 50 M0?/r~\ 
their luminous counterparts {M{z) < —22) shows SFRs greater than 300 and as high 
as 2000 MQyr~^ , suggesting that luminous and ultra- luminous infrared galaxies in the 

■ low-redshift Universe could be the progenitors of massive nearby E+A galaxies. We 
\ perform a comprehensive study of the characteristics of the quenching that truncates 
i the starburst in E+A systems. We find that, for galaxies less massive than IO^^Mq, the 

CO ' quenching efficiency decreases as the galaxy mass increases. However, for galaxies more 

massive than IO^^Mq, this trend is reversed and the quenching efficiency increases with 

■ galaxy mass. Noting that the mass threshold at which this reversal occurs is in excellent 
' agreement with the mass above which AGN become significantly more abundant in 
, nearby galaxies, we use simple energetic arguments to show that the bimodal behaviour 

of the quenching efficiency is consistent with AGN and supernovae (SN) being the 
principal sources of negative feedback above and below M ~ IQ^^Mq respectively. 
^ . The arguments assume that quenching occurs through the mechanical ejection or 

\—{ ' dispersal of the gas reservoir and that, in the high mass regime (M > IO^^Mq), the 

5^ \ Eddington ratios in this sample of galaxies scale as M'^ , where 1 < 7 < 3. Finally, we 

use our E+A sample to estimate the time it takes for galaxies to migrate from the 
blue cloud to the red sequence. We find migration times between 1 and 5 Gyrs, with 
a median value of 1.5 Gyrs. 

Key v^rords: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: fundamental parameters 
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1 INTRODUCTION 

'E+A' galaxies exhibit strong Balmer absorption lines, char- 
acteristic of significant recent star formation, but lack [Oil] 
or Ha emission which are characteristic of ongoing star 
formation. Their spectral features indicate that a vigor- 
ous episode of recent star formation in these systems has 
been quenched abruptly. As 'post-starburst' systems, E+A 
galaxies provide a valuable evolutionary link between the 
gas-rich star-forming population and gas-poor quiescent sys- 



tems. The persistent dichotomy in galaxy colours over a 
large range in redshift, coupled with the gradual build-up of 
the red sequence over time (e.g. iBell et al]|2004 ). indicates 
that how and over what timescale galaxies evolve from the 
'blue cloud' to the 'red sequence' is of significant importance 
in understanding the macrosopic evolution of the galaxy 
population over the lifetime of the Universe. The charac- 
teristics of E+A galaxies provide a unique (but transient!) 
window into the very point where this transition from blue- 
cloud to red sequence is about to begin. Understanding what 
drives their vigorous star formation episode and, in particu- 
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lar, what causes it to be quenched so suddenly is clearly an 
important step in o ur understanding of gala xy evolution. 

First noticed by Dressi er fc Gunn * (1983) in distant {z ~ 
0.4) clusters, E+A galaxies have been subsequently found to 
be abundant in all types of environments. At intermediate 
redshifts (0.3 < z < 1), the proportion of E+A galaxies ap- 
pears to_be_a^factorof four greater in clusters than in the 
field (|Tran et al.l l2004'l . In the nearby Univ erse E+A galax- 
ies h a ve been detected ma inly in the field l|Zabludoff et al.l 
1 19961 : iQuintero et al.l 120041 ') primarily because most galax- 
ies do not reside in dense cluster-like environments. Indeed, 
the overall distribution of the local environments of E+A 
systems follows that of the general galaxy population as a 
whole (Blake et al. 2004). The fraction of E+A systems in 
clusters shows a rapid decline from intermediate redshifts 
(z ^ 0.5) where it is ty pically higher than 20 percent (e.g. 
ICouch fc Sharples|[l987 : Bclloni ct al. 19 9^) t o less than 1 
percent in local clusters (jFabricant et anil99ll ). 

Several plausible explanations for the E+A phe- 
nomenon have been proposed and studies indicate that there 
are multiple channels for creating such post-starburst spec- 
tral signatures. Since E+A galaxies were first detected in 
clusters, it was initially thought that their production re- 
quired a cluster-specific mechanism, such as galaxy harass- 
ment or ram-pressure stripping. However, their abundant 
presence in the field indicates that o ther channels exist for 
the production of E+A systems (e.g. iGot a 2005 1. although 
cluster-specific mechanisms may contribute or even dom- 
inate their evolutionary pathways in dense regions of the 
Universe. 

Another possible explanation for E+A-like spectra is 
that the optical e mission lines are heavily suppressed by dust 
obscuration (e.g. lSmail et al.ll 19991 : ICouch fc Sharples|[l987l : 
iPoggianti fc Wull2000l ). An efficient test of this scenario is 
to observe the E+A population in the radio wavelengths 
where the effects of dust obscuratio n are absent. However , 
a radio st udy of 15 galaxi e s in t he IZabludoff et al.l (|l996l ) 
sample by iMiller fc OwenI (120011 ) detected only two at ra- 
dio luminosities which are consi s tent w ith modest levels of 
star formation. Similarly, iGotd (120041 ') performed VLA 20 
cm radio-continuum observations of 36 E+A galaxies drawn 
from the SDSS and found that none of these galaxies were 
detected in the radio wavelengths. These results are sup- 
ported by'Galaz' (2000) who do not find compelling evidence 
for strong internal dust extinction in their E+A sample. 

Some E+A galaxies exhibit morphological d isturbances 
indica tive of interactions with near neighbours. lYang et al.l 
l|2004f ) have studied the m orphologies of th e five b luest E+A 
galaxies in the sample of IZabludoff et al] (|l996l ) using the 
WFPC2 camera on board the Hubble Space Telescope (HST) 
and found that four out of these five galaxies display mor- 
phological disturbances consistent with recent mergers. Sim- 
ilarly, Goto (2004) have shown, using the largest sample of 
E+A galaxies studied to date, that these systems have an 
excess of local galaxy density at spatial scales less than 100 
kpc but not at scales spanned by galaxy clusters (~ 1 Mpc) 
or large-scale structure (~ 8 Mpc). They find that ~ 30 per- 
cent of E+A galaxies exhibit morphological disturbances, in- 
dicating that their evolution is linked, at least partially, to 
me rgers and in teractions. It is worth noting that, although 
the lGotol (|2004l ') sample is larger and their results are statisti- 
cally more robust, detection of small, low-surface brightness 



(tidal) features probably requires both deeper and higher 
resolution imaging than that provided by the SDSS. Hence, 
the fraction of E+A galaxies w ith tid al tails may well be 
higher than ~ 30 percent in the lGotol (,2004 ) sample, mak- 
ing the conclusions from these two studies more consistent. 

From a theoretical perspective, the merger hypothesis 
is supported by numerical simulations which indicate that 
gas-rich mergers are capable of triggering strong star for- 
mation episodes. The immediate aftermath of such an event 
is predicted to be a largely spheroidal remnant with a blue 
core, widespread morphological disturbances such as tidal 
tails and young st ar clusters (e.g. Mihos fc Hcrnquist 199^; 
iBekki et al]|200ll '). Indeed, the study of IVang et al.l (|2004h 
has found a striking correspondence between the predicted 
properties of merger remnants and the observed character- 
istics of blue E+A galaxies. 

Numerous authors have studied E+A galaxies using a 
variety of indicators, which allows us to infer the charac- 
teristics of the (recent) starburst that dominates the spec- 
tra of these systems. 'While the age of the burst should be 
approximately ~ 1 Gyr, simply by virtue of A-type stars 
being present, the mass f ractions forming in the burst are 
inferred to be quite high. lYang et al.l (|2004h estimate that 
50 percent of the mass may form in the burst, allowing the 
E+A galaxies to fade on to the E/SO fundamen tal plane in 
a few Gyrs aft er the event. Norton et al.l ([iooj), who anal- 
ysed the entire IZabludoff et al.l ( 19961 ) sample, estimate that 
the fraction of stellar mass in yo ung (A-type) stars is typi- 
cally between 30 and 80 percent. iBressan et all (|200ll ) infer 
a similarly high mass fraction of 60 perc ent for the sample of 
E+A galaxies in IPoggianti et all (|200ll ). while iLiu fc GreenI 
(j 19961 ) derive mass fractions typically greater than 20 per- 
cent (and as high as 100 percent) for eig ht E+A galaxies 
taken from the literature (but note that iLiu et alT (|2007l ) 
revise the inferred mass fraction to lower values for one of 
their E+A galaxies - G515). 'While E+A's are, by definition, 
presently quiescent, an estimate of the star formation rate 
(SFR) that might have characterised the recent starburst 
can be estimated from dusty starburst galaxies which show 
an E+A signature in their spectra. .Poggianti fc "Wu ( 200(3) 
estimate the SFRs in these systems to be between 50 and 
300 M0 per year. 

'While mergers and interactions might trigger star for- 
mation, the exact mechanism by which the starburst is sub- 
sequently quenched remains largely unexplained. Star for- 
mation is, of course, inherently subject to 'negative feed- 
back', since the star formation rate (SFR) depends on the 
mass of cold gas available. As star formation proceeds and 
the gas reservoir is exhausted, the rate of further star forma- 
tion activity correspondingly decreases. However, the prop- 
erties of E+A galaxies are generally inconsistent with such 
a gradual, self-modulated decline in the SFR and imply a 
more abrupt quenching of star formation activity. 

Common sources of kinetic and thermal 'feedback' that 
might be capable of quenching star formation, e.g. by eject- 
ing or heating the gas reservoir, include supernovae (SN), 
which are effective mainly in low-mass galaxies due to their 
shallow potential wells, and AGN, which are thought to op- 
erate in massive galaxies where SN feedba ck is not ene rgetic 
enough to 'disable' the cold gas reservoir. iGotd (|2006D have 
demonstrated that some galaxies which contain signatures 
of an active AGN also exhibit post-starburst signatures, in- 
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Figure 1. The position of our E+A galaxies (filled blue circles) 
in the (NUV — r) vs (g-r) colour space, compared to a sample of 
early-type galaxies drawn from the SDSS DR5 using the method 
of Kaviraj et al. (2006b) and cross-matched with GALEX. The 
E+A population inhabits the blue cloud and is well separated 
from the early-type galaxies which define the bulk of the red se- 
quence. Note that the colours are in the observed frame i.e. they 
are shown without applying K-corrections. 




2=3 

(ti ~ 10 Gyrs) 

Figure 2. Model SFHs (see Section 2) are constructed by assum- 
ing that an instantaneous burst of star formation at high redshift 
{z ~ 3) is followed by a second burst, modelled as an exponen- 
tial. The free parameters are the age (t2), mass fraction (/2) and 
timescale (ra) of the second burst. t2 is allowed to vary from 
0.05 Gyrs to the look-back time corresponding to ^ = 3 in the 
rest-frame of each E+A galaxy. /2 varies between and 1 and T2 
is allowed to vary from 0.01 Gyrs to 4 Gyrs. r^yn indicates the 
dynamical timescale of the galaxy. More catastrophic quenching 
leads to smaller values of the ratio T2 / T^y^^ , which is used in Sec- 
tion 5 to parametrise the efficiency with which quenching occurs 
in the E+A galaxies studied here. 



dicating that a connection exists between these two phenom- 
ena. They find that the post-starburst regions are centred 
around the AGN, the spatial proximity implying that AGN 
activity could plausibly affect the region around it. Indeed, 
the fact that the AGN outlives the starburst is an indica- 
tion that it may have played a role in the quenching of star 
formation. 

In this study we extend the existing literature on post- 
starburst galaxies by performing the first study of E+A sys- 
tems that incorporates their ultraviolet {UV) photometry. 
The extreme sensitivity of the UV to young stars allows us to 
put strong constraints on the star formation history (SFH) of 
galaxies within the last ~ 2 Gyrs. The addition of UV pho- 
tometry alleviates the age-me tallicity degeneracy that com- 
monly plagues optical studies l|Kavirai et al.|[2"007^ . allowing 
us to better quantify the recent SFH than can be achieved 
using optical data alone. Indeed, the sensitivity of the UV to 
young stars has recently been been exploited to demonstrate 
that early-type galaxies, traditionally thought to be entirely 
quiescent systems, exhibit widespread evidence for low-level 
recent star formation, consistent with the expectations of 
the standard LCDM model (Kaviraj et al. 2006a). 

In this paper, we exploit the properties of the UV to 
accurately reconstruct the recent star formation history of 
individual E+A galaxies. By comparison to a large library of 
models we derive estimates for the ages of the recent bursts, 
the mass fractions formed in them and the timescales over 
which they occurred. By comparing the derived timescales 
to the dynamical timescales of the galaxies, we explore the 
efficiency with which star formation has been quenched in 
these systems and study how this quenching efficiency varies 
as a function of galaxy properties such as mass, luminosity 
and stellar velocity dispersion. We use simple energetic argu- 
ments to study whether the characteristics of the quenching 
are consistent with feedback from supernovae and AGN. Fi- 



nally, we study the timescale over which E+A galaxies can 
evolve from their present positions in the blue cloud to the 
red sequence. Since E+As are, by definition, at the point 
of truncation of their star formation, it is possible to get a 
rather 'clean' estimate of this migration timescale by study- 
ing the evolution in both UV and optical colours. 



2 SAMPLE SELECTION 

We use the publicly available Garching SDSS catalog to se- 
lect E+A galaxies from the S DSS DR4. T he criteria used for 
E+A selection follow that of lGotol (|2005l ): H(5 (EW) > 5.0A 
(measured in the wider window of 4082-4122A), Ha (EW) 
> -3.0A, [OH] (EW) > -2.5l. Absorption lines have a 
positive sign in these definitions. Note that including the 
equivalent width of the Ha line makes the selection criteria 
stronger, because using only the [OH] line could result in 
more than 50 percent of the selected galaxies having measur- 
able Ha emission (e.g. Goto et al. 2003). In addition to the 
line measurements we also use estimates of stellar mass and 
AGN classifications provided by the Garching SDSS catalog. 
The line index m easurements in this cat alog were calculated 
using the code of iTremonti et al] (|2004l ). while stellar mass 
estimates and AGN classifications were taken from a series 
of comprehensive analyses of the local galaxy population ob- 
served by the SDSS (Kauffmann et al. 2003a; Kauffmann et 
al. 2003b; Brinchmann et al. 2004) 

To ensure the accuracy of the line measurements, the 
sample was restricted to galaxies where the median S/N 
of the spectrum is greater than 10. Galaxies which show 
evidence of an active AGN were removed, because scattered 
light from the AGN could contaminate the UV continuum 
and affect the derived parameters presented later in this 

^ http:/ /www. mpa-garching.mpg.de/SDSS/DR4/ 
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analysis. Finally, the sample was restricted to galaxies in 
the redshift range < z < 0.2, since the UV filters used 
in the analysis (see below) trace flux blueward of the UV 
continuum at higher redshifts. 

This E+A sample, drawn from the SDSS DR4, was 
then cross-matched with publicly available UV photome- 
try from the second data release of the GALEX mission 
(Martin et al. 2005). GALEX provides two UV filters: 
the far- ultraviolet (FUV), centred at ~ 1530A and the 
near-ultraviolet (NUV), centred at ~ 2310A. This cross- 
matching produced 38 E-l-A galaxies which have at least a 
detection in the NUV filter, 28 of which have photometry 
in both the FUV and NUV filters. 

Figure [1] shows the position of our E+A galaxies (filled 
blue circles) in the {NUV ~ r) vs {g — r) colour space, com- 
pared to a sample of early-type galaxies drawn from the 
SDSS DR5 using the method of Kaviraj et al. (20Q6b) and 
cross-matched with the GALEX second data release. It is 
apparent that the E-l-A population is well separated from 
the early-type population which defines the bulk of the red 
sequence. Comparison to Figure 1 in Yi et al. (2005) indi- 
cates that E-l-A galaxies inhabit the blue cloud. 

Since previous studies have indicated that E-f A galax- 
ies are potential progenitors of early-type galaxies, we check 
the morphologies of the E-l-A galaxies in our sample using 
the SDSS fracDev parameter. The SDSS pipeline fits both 
deVaucouleur's and exponential surface brightness profiles 
to galaxy images (in each filter) and creates a composite 
best-fit profile using a linear combination of the two fits. 
The fracDev parameter is the weight of the deVaucouleur's 
profile in this composite fit. Thus, galaxies with large val- 
ues of fracDev can be considered to have spheroidal mor- 
phology. We find that all but 6 of the galaxies studied in 
this sample have fracDev > 0.8 in the r-band, which in- 
dicates that the majority of E-|-A galaxies in this sample 
are indeed spheroidal systems. It is worth noting that all 6 
galaxies which have fracDev < 0.8 are low-mass systems 
(M < 10^° Mq) and their non-deVaucouleur's morphologies 
are consistent with the fact that dwarf ellipticals typically 
show exponential surface brightness profiles. 



3 PARAMETER ESTIMATION 

We estimate parameters governing the star formation his- 
tory (SFH) of each E-fA galaxy by comparing their 
{FUV, NUV, u,g,r, i, z) photometry to a library of synthetic 
photometry generated using a large collection of model 
SFHs. Each SFH is constructed by assuming that an in- 
stantaneous burst of star formation at high redshift {z ~ 3) 
is followed by a second burst which is modelled as an expo- 
nential. 

Figure [2] shows a schematic representation of the model 
SFHs. The free parameters in this analysis are the age {t2), 
mass fraction (/2) and timescale (T2) of the second burst. t2 
is allowed to vary from 0.05 Gyrs to the look-back time cor- 
responding to z = 3 in the rest-frame of each E-l-A galaxy. 
/2 varies between and 1 and T2 is allowed to vary from 
0.01 Gyrs to 4 Gyrs. Since our focus is purely on E-|-A 
galaxies, we must, by definition, exclude models which con- 
tain ongoing star formation. Therefore, we only keep mod- 
els where the intensity of star formation in the second burst 



at present-day is less than 5 percent of the intensity when 
the burst startetfl. Changing this threshold to less than 5 
percent leaves our results unchanged. Note that excluding 
models that have ongoing star formation reduces the allowed 
parameter space of model SFHs, resulting in much tighter 
constraints on the parameters, in particular, the timescale 
of star formation T2. 

The motivation for an instantaneous burst at high red- 
shift stems from the fact that, as potential progenitors of 
spheroidal galaxies, one expects a substantial fraction of 
stell ar mass to have for med at high redshift in these galaxies 
fe.g. iBower et al.lll992l . Kaviraj et al. 2006b). Furthermore, 
since the UV is insensitive to stellar populations older than 
~ 2 Gyrs coupled with the fact that optical c olour evolution 
virtually stops after 5-6 Gyrs (e.g. IyI 120031 ). resolving this 
old burst does not affect the UV and optical colours of the 
model. 

The UV flux in each model is completely dominated by 
the second burst of star formation. Even though our expec- 
tations, on the basis of the properties of E-l-A galaxies (and 
previous studies of E-l-A systems) , are that t2 ~ 1 Gyr and 
/2 is high, our model library does not put a prior on these pa- 
rameters - vigorous starbursts at recent epochs are allowed, 
as is slowly declining star formation from intermediate red- 
shifts. In essence, the comparison between the synthetic li- 
brary and observed photometry is allowed to 'choose' the 
optimum SFH that satisfles the photometry of each E-l-A 
galaxy. 

To build the library of synthetic photometry, each 
model SFH is combined with a single metallicity in the range 
O.IZ0 to 2.5Zq and a value of dust extinction parametrised 
by E(B-V) in the range to 0.5. Photometric predictions 
are generated by combining each model SFH with the cho- 
sen metallicity and E(B- V) values and convolving with the 
stellar models of^ (1200^ ) through the GALEX FUV, NUV 
and SDSS u, g, r, i, z fllters. This procedure yields a synthetic 
library of 1.5 million models. 

Since our E+A galaxies are observed at a range of red- 
shifts, equivalent libraries are constructed at redshift inter- 
vals of (5z = 0.01. A fine redshift grid is essential in such a low 
redshift study because a small change in redshift produces 
a relatively large change in look-back time over which the 
UV flux can change substantially, inducing 'K-correction- 
like' errors into the analysis. 

The free parameters {t2, f2 and T2) are estimated by 
comparing each observed galaxy to every model in the syn- 
thetic library, with the likelihood of each model (exp — x^/2) 
calculated using the value of computed in the standard 
way. From the joint probability distribution, each param- 
eter is marginalised to extract its one-dimensional proba- 
bility density function (PDF). We take the median of this 
PDF as the best estimate of the parameter in question and 
the 16 and 84 percentile values as the 'one-sigma' uncertain- 
ties on this estimate. The cosmological parameters used in 
this study assume a ACDM model: h = 0.7, = 0.3 and 
= 0.7. 

Note that there are 10 galaxies in our sample which 
do not have FUV detections. The parameter estimation on 



^ This implies that models with small values of 12/ T2 are ex- 
cluded. 
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Figure 3. TOP PANEL: Ages and mass fractions of the burst 
ttiat creates tlie E+A galaxies in our sample. Burst ages are typ- 
ically within a Gyr, as one expects from the lifetimes of A-type 
stars which contribute significantly to E+A spectra. We recover 
large mass fractions, typically higher than 10 percent and possibly 
as high as 70 percent, consistent with previous studies. BOTTOM 
PANEL: Effective timescales of the starbursts in E+A galaxies. 
The timescales are typically short - between 0.01 and 0.2 Gyrs - 
implying that the star formation rates during the burst phase are 
high (see text in Section 4 and Figure |4]l. 



these galaxies is carried out using only the NUV, u,g,r,i 
and z-band filters. We do not use the FUV detection limit 
as an upper limit to the FUV flux. 



4 CHARACTERISTICS OF THE RECENT 
BURST: AGES, MASS FRACTIONS AND 
TIMESCALES 

We begin by presenting estimates for the basic parameters 
(fa, /2 and T2) that describe the burst of star formation that 
immediately precedes the quenching in our E+A galaxies. 
The top panel of Figure [3] shows the ages and mass fractions 
of the recent burst in each E+A galaxy. The burst ages are 
typically within 2 Gyrs, as one expects from the lifetimes of 
A-type stars which contribute significantly to E+A spectra. 
In agreement with previous studies, we derive large mass 



Figure 4. Implied SFRs during the burst phase in E+A galaxies. 
The implied SFR is estimated by dividing the stellar mass formed 
by the timescale of the burst. Note that, since the burst lasts 
longer than one timescale, the SFRs are overestimated. However, 
the estimates probably give a reasonable indication of the SFR 
during the most intense period of the starburst. 



fractions, which are typically higher than 10 percent and 
possibly as high as 70 percent in some galaxies. The bottom 
panel shows the timescales over which the bursts take place. 
The timescales are typically short - between 0.01 and 0.2 
Gyrs - implying that the star formation rates during the 
burst phase are extremely high. 

We can estimate this 'implied' SFR by dividing the stel- 
lar mass formed by the timescale of the burst (Figure |4}. We 
note that, since the burst leists longer than one timescale, 
the SFRs are overestimated. However, the estimates prob- 
ably give a reasonable indication of the SFR during the 
most intense period of the starburst and we find a tight 
correlation between the mass of the E+A galaxy and the 
implied SFR. (see bottom panel of Figure |3} . While low- 
lumnosity {M{z) > -20) E+As have implied SFRs of less 
than 50 M^yr'^ , E+A systems at the high-luminosity end 
{M{z) < -22) exhibit SFRs of greater than 300, and as high 
as 2000 MQyr'^. 

The high SFRs inevitably lead to comparisons 
with luminous and ultra-luminous infrared galax- 
ies (LIRGs/ULIRGs), whose high infrared luminosi- 
ties {LiR > lO^'^Lp)) imply mas sive ongoing star- 



bursts (e.g ISanders fc Mirabell 1 19961 ). SFRs in LIRGs 



(lO^^Z/Q < LiR < IQ^^Lq) typically exceed a few tens of 
solar masses per year, while in ULIRGs {Lir > IO^^Lq) 
the SFRs can be as high as hundreds of solar masses 
per year (e.g. lKennicutt) [l998l ). IWang et all (|2006D . who 
have performed a study of LIRGs in the nearby Universe 
{z ~ 0.1), find SFRs of between 10 and 100 M^yr'^ for 
their sample of galaxies (see their Figure 5). Noting that 
our SFRs are overestimated, and could be a few factors too 
high, we suggest that LIRGs at low redshift could transform 
into massive E+A galaxies. However, we note that the 
robustness of such a conclusion depends on performing an 
identical parameter estimation on LIRGs at low redshift 
and by considering whether other factors, such as the local 
environments of E+A galaxies, correspond closely to those 
of LIRGs. 
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Figure 5. The timescale ratio, = T2/T^y^, plotted against 
galaxy mass (top), absolute z-band luminosity (middle) and stel- 
lar velocity dispersion (bottom). Note that the red dotted lines 
indicate standard least-squares fits to the data. The fits are per- 
formed separately for the massive (M > IQ^^Mq) and low -mass 
(M < lO^^Mo) E-l-A galaxies in each panel. 



5 THE EFFICIENCY OF QUENCHING AS A 
FUNCTION OF GALAXY MASS 

One of the primary goals of this paper is to study the char- 
acteristics of the quenching that truncates the burst and 
determine the possible mechanisms by which this quenching 
could occur. A key parameter in this analysis is the efficiency 
with which the star formation is quenched. We study this 
'quenching efficiency' by comparing the effective timescale 
of the burst, T2, to the dynamical timescal^fl of the galaxy, 
Tdyn, which is defined as 



\GM) ' 



(1) 



and describes the 'natural' timescale over which pro- 
cesses such as star formation would take place if left un- 
hindered 0. Note that, in the calculations that follow, we 
have used the stellar mass of the galaxy as the value of M. 
We parametrise the quenching efficiency by constructing the 
timescale ratio 



Tt = T2/Tdyn- 



(2) 



More efficient quenching reduces T2 and therefore pro- 
duces a lower value for Tt- The timescale ratio is, therefore, 
inversely correlated with quenching efficiency - more catas- 
trophic quenching results in lower values of r^. 

In Figure [5] we plot against three independent indi- 
cators of galaxy mass/luminosity. The galaxies are colour- 
coded using their {NUV — r) colour. The plot of r-r vs galaxy 
mass (top panel), reveals a striking difference in the be- 
haviour of rr with increasing galaxy mass. We find that, be- 
low a mass of IO'^'^Mq, the quenching efficiency decreases 
with increasing galaxy mass. However, when the galaxy 
mass is greater than ~ lO^^M© this trend is reversed and 
the quenching efficiency then correlates positively with the 
galaxy mass. We check this behaviour with both M{z) and 
stellar velocity dispersion. We find that the correlation of 
r-r with galaxy mass and, in particular, its apparent re- 
versal is present in all the separate indicators of galaxy 
mass/luminosity. While a strong trend is apparent with stel- 
lar mass (M) and z band luminosity, the trend is weaker in 
the stellar velocity dispersion (a), most probably driven by 
bigger uncertainties in the a measurement (as suggested by 
the larger scatter in this plot). A plausible source of un- 
certainty is the fact that the E-l-A galaxies are unrelaxed 
(indeed a significant fraction look disturbed in the SDSS 
images), implying that a does not yet trace the virial mass 
of the galaxy. An additional source of uncertainty is the pres- 
ence of a substantial young stellar population. The presence 
of young sub-components could create a bias towards lower 
values of o, with the magnitude of the bias driven by the 
mass fraction and age of the young stars and on the exact 
size of the physical aperture subtended by the SDSS spectra 
in each galaxy. We therefore suggest that the correlation be- 



^ Using s = ^at^ from the equations of motion, the dynamical 
timescale follows by setting a = g = (GM/R?) and s = R. 
^ Note that we use the radius that contains 90 percent of the 
Petrosian flux in the i-band as the value of R. 
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tween r-r and a should be treated with caution considering 
the caveats we have presented above. 

Figure [5] shows that the reversal in the correla- 
tion occurs at M ~ 10^°Mq, M{z) ~ -21 and a ~ 100 
kms~^ when Tt is plotted against mass, absolute z-band lu- 
minosity and stellar velocity dispersion respectively. The di- 
chotomy in the trend of r-r against galaxy mass indicates 
that the principal quenching mechanisms for galaxies above 
and below the mass threshold of IO^^Mq behave very dif- 
ferently. In particular, the mechanism that operates in the 
regime M < 10^'' Mq becomes weaker as the galaxy mass 
increases while the mechanism that operates in the regime 
M > 10^° AIq becomes stronger as the galaxy mass in- 
creases. 

The mass threshold (10^"Mq) at which the vs. M 
trend reverses is particularly interesting because it is in ex- 
cellent agreement with the mass above which AGN become 
significantly more abundant in galaxies in the nearby Uni- 
verse (Kauffmann et al. 2003a). Given this result, the re- 
versal at M ~ 10^° Mq can be explained, at least qualita- 
tively, in terms of SN and AGN being the principal quench- 
ing sources above and below this mass threshold. In the ab- 
sence of AGN, the primary source of negative feedback are 
SN. As galaxies become more massive and the depth of the 
potential well increases, SN find it increasingly more difficult 
to eject gas from the system. This results in the quenching 
becoming less efficient at higher galaxy masses, exactly as 
we observe in Figure (5] 

However, once AGN begin to appear (above M ~ 
10^" Mq), they become the dominant source of negative feed- 
back. For a 10^^ Mq galaxy, the energy input from a popu- 
lation of (Type II) SN formed as a result of a starburst that 
creates half the stellar mass of the galaxy is ~ 3 x 10^^ erg^, 
while the energy output from a central AGN, emitting at the 
Eddington limil[f| over a timescale of 0.1 Gyrs is ~ 3 x 10^* 
ergfl In terms of energy input the AGN dominates and it 
is reasonable to assume that it dictates the mechanics of the 
quenching. Note that we did not consider Type la super- 
novae in this calculation because the derived timescales (see 
Figure [3]) are considerably shorter than 1 Gyr, which is the 
typical time delay between the onset of star formation and 
the appearance of Type la supernovae. 

Given that the mass of the black hole (Mbh) scales 
with the centra l velocity dispersion a as Mbh ~ o"'' 
IIGebhardt et al ] |200Q : iFerrarese fc Merrittii2000h , where P 



varies in range 4-5, and assuming the energy outputted by 
the AGN scales with Mbh (c.f. the Eddington luminosity 
Ledd oc Mbh), we expect AGN feedback to become more 
effective, and thus the quenching efficiency to increase, as 
galaxy mass increases. Again, this is indeed what we ob- 
serve in Figure O We briefly note that, by definition, our 
sample does not contain galaxies that show signs of current 
AGN activity, implying that the AGN has shut itself off in 
the process of providing feedback. If, as we argue below, me- 
chanical feedback from the AGN depletes the gas reservoir, 
then this process also removes the AGN's own fuel source. 
Thus it is plausible that the feedback process simultaneously 
quenches both star formation and AGN activity. 

It is worth noting that a similar dichotomy in 
the primary feedba ck mechanism has been inferred by 
l|Shankaret al.ll2006l ). from a study of the stellar and bary- 
onic mass functions of galaxies, extracted using the mass-to- 
light ratios of stars and gas derived from galaxy kinematics. 



6 THE EXPECTED DEPENDENCE OF THE 
QUENCHING EFFICIENCY ON GALAXY 
MASS 

We now use simple energetic arguments to derive the ex- 
pected relationship between the timescale ratio r-r and 
galaxy mass in the two cases where the quenching is due 
to SN and AGN respectively. 

We make the assumption that quenching takes place 
due to the ejection of available gas from the potential well. 
The energy required for quenching, Eq, is given by, 



P _ 1 2 



(3) 



where rrig is the mass of gas to be ejected and v^^^ is 
the escape velocity. Since, 



this implies that 



Eq oc rrig. — , 



(4) 



(5) 



or equivalently, in terms of parameters derived in this 
analysis that. 



^ The energy injected by a population of SN is given by 
VSNssn^M , where r?s]v(~ 5 X 10~^Mq"') is the number of SN 
per unit stellar mass and esjv(~ 10^^ ergs) is the energy supplied 
by each supernova. The fraction of that energy that couples me- 
chanically to the material in the ISM is approximately given by 
(t/20, 000 kms~^, where a is the velocity dispersion and 20,000 
kms~^ is the velocity of the SN ejecta. For a M 10^^ Mq 
galaxy, a ~ 200 kms~^. 

^ The Eddington luminosity for a mass M is ~ 1.3 X IO'^^M/Mq 
ergs per second. 

A A/ ~ lO^^Af0 typically hosts a M ~ lO^M p) black hole, 
since (MBH/Mtuige) ~ 10"^ tearing fc Rixll2004) . The Edding- 
ton luminosity of such a black hole is ~ 10*® ergs , of which a 
fraction cr/c kms~^ couples to the mechanically to the material 
in the ISM. a is the velocity dispersion of the galaxy and c is the 
speed of light. 



/ M \2/3 /Mr^\2/3 

Eq(Xmg.[ ) oc TTlg. ) . (6) 

^Tdyn' \ T2 / 

The relationships in Eqn (6) use substitutions from 
Eqns (1) and (2). Eqn (6) describes the dependence of the 
energy required for quenching on the timescale ratio r-r, the 
burst timescale T2 and galaxy mass M. 

We now make the simplifying assumption that the mass 
of gas (nig) required to be ejected for quenching is approx- 
imately constant in all galaxie s. This is con sistent with the 
downsizing phenomenon (e.g ICowie et al.l [l996. ') , whereby 
smaller galaxies exhibit higher star formation rates than 
their massive counterparts at lower redshifts. Given the ap- 
parent universality of the Schmidt-Ken nicutt law for (qui- 
escent) star formation (|Kennicuttlll998l ). downsizing implies 
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that smaller galaxies tend to be more gas-rich than larger 
ones, justifying the empirical assumption that their abso- 
lute neutral gas reservoirs may have similar siz es. Cold gas 
fract ions in low-redshift LIRGs (see Figure 5 in I Wang et alj 
I2OO6I ). which are potential progenitors of E-|- A galaxies (Sec- 
tion 4), indeed show a declining trend with galaxy mass. 
Across the mass range 

IQiO-S _ IO^'^Mq, the gas fractions 
increase from less than 5 percent to greater than 25 percent, 
suggesting that the cold gas reservoirs in these galaxies are 
similar in size and that our assumption of constant mg is 
a reasonable one. A small spread in the real values of 
will induce a scatter in the trends derived in Section 6.1 and 
6.2 below, without perturbing the correlations between the 
variables. With this assumption Eqn (6) becomes 



Eq OC I 



2/3 ^ 



(7) 



We should also note that the assumption that quench- 
ing occurs through the mechanical ejection of gas may im- 
ply that E-f A galaxies should be relatively gas-poor, al- 
though the gas could simply be dispersed away from the 
central regions of the potential well where star formation 
would normally take place. Two studies of th e HI content 



of_E+A galaxies have been conducted so f ar. IChang et al 



2001 ) observed five E-l-A galaxies from the lZabludoff et"al 



19961 ) sample and detect ed only one with a n HI mass of 
~ 3.4 X 10^ M0. Similarlv. iBuvle et all (120061 ) detected sim- 
ilar amounts (> 10^ M0) of HI in three out of six E-|-A 
galaxies they observed in HI. While , it might be us eful to 
compare the gas masses estimated bv lBuvle et al.l fe.g |2006l ) 
for E-l-A galaxi es to those es t imate d for LIRGs (as potential 
progenitors) by IWang" et al.1 (|2006l ). the mass ranges traced 
by these two studies are quite different. iBuvle et af ('2006|) 
trace galaxies with M < lO^°-^Me0, while Wang c t'aLl 
120061 ') study galaxies with M > 10^° '^ Mq). The small over- 
lap region around 10^"''' Mq indicates that the HI fractior(f| 
in the LIRGs may be higher by a factor of 3. We note, how- 
ever, that due to the inadequate overlap between the two 
samples, this result is not very robust. However, given the 
apparent lac k of H I in the majority of E-|-A's observed by 
Chang et al.l (I2OOII ) and the smaller gas fractions found in 
Buvle et al.l l|2006l ) compared to nearby LIRGs (which have 
comparable SFRs to our E-l-A galaxies) the assumption of 
quenching through the mechanical ejection of gas seems rea- 
sonable. 
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Figure 6. vs. Qsn for galaxies with M > IO^^Mq (black) 
and A/ < IO^^A/q (red). The solid lines indicate standard least- 
squares fits to the data. 



E^^ = 77siv£SivAM, 



(8) 



where r]sN is the number of SN per unit stellar meiss, 
esN is the energy supplied by each supernova and AM is 
the amount of stellar mass formed. Hence, 



E^^ OC /2M, 



(9) 



since AM = /2M. Combining Eqn (7) and Eqn (9) 



gives. 



so that 



log(r|/|M) 



21og(rr) + k 



(10) 



(11) 



As we would naively expect from SN driven feedback, 
Tt correlates with the galaxy mass (M). Since is inversely 
correlated to the quenching efflciency, this implies that, as 
galaxy mass increases, the quenching efficiency decreases in 
the SN feedback scenario. Defining Qsn = r|/|M we have. 



6.1 Quenching by supernovae 

The energy injected by SN into the ISM is proportional to 
the amount of stellar mass produced in the burst so that 



* The masses of the lBuvle et al.l ll2006t ) galaxies have been esti- 
mated by converting the values of M{B) given in their Table 1 to 
M{z), which are then converted to masses using the M vs. M(z) 
relation for the E-l-A galaxies in our sample. 

^ Note that we h ave converted the H 2 mass (which is the quan- 
tity estimated bv lWane et al.l | |200(J) ) to an HI mass, using the 
average HI/H 2 ratio for late-type galaxies (~ 0.6) given by 
iFukugita et al.l lll998l . see their Table 1) 



log(Qsiv) = 21og(rO + A;, (12) 

where k absorbs the constants in the proportionality 
relations used to derive Eqn (12). A plot of log(Qs]v) against 
log(rT-) should therefore show a gradient of 2, if the systems 
being considered are indeed quenched purely by SN. 

In Figure |S] we plot vs. Qsn for both low-mass 
(M < 1O^°M0; black) and massive (A/ > 1O^°M0; red) 
galaxies. The solid lines indicate standard least-squares fits 
to the data. We find that, while low-mass galaxies show a 
gradient of 1.98 ±0.18 (and are therefore consistent with the 
expected gradient of 2), the best-fit relation for the massive 
galaxies is too shallow (0.56 ±0.06) and does not satisfy the 
predictions. We conclude, therefore, that low-mass galaxies 
show good consistency with being quenched purely by SN 
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Migration time to red sequence (Gyrs) 

Figure 7. Histogram sliowing tlie migration time of E+A galaxies 
from the blue cloud to the red sequence in both the UV and 
optical colours. The median migration time (indicated by the solid 
red lino is ~ 1.5 Gyrs. The y-axis shows the frequency normalised 
to 1.) 



feedback while massive galaxies are essentially inconsistent 
with such a picture. 



6.2 Quenching by AGN 

The energy injected by an AGN (E^'^^) is driven by its 
luminosity (L) so that, 



TpAGN r 

bj OC L.T2. 



(13) 



Rewriting Eqn (13) in terms of the Eddington luminos- 
ity (Ledd), we have 



■T2, 



(14) 



where p = L/Ledd is the Eddington ratio. L^dd is pro- 
portional to the mass of the central black hole {Mbh), which 
in turn has been shown to scale with the central velocity 
dispersion (cr) of the galaxy as Mbh = ■ E stimates of 
P vary in the range 4-5, with Gebhardt et al.l (e.g . 
reporting j3 = 3.75 ± 0.3, while iFerrarese fc MerrittI 
report (3 — 4.8 ± 0.5. For simplicity, we begin by assuming 
that ^ is constant across our entire sample of AGN. Then, 



2000ll 



20001 ) 



where we have used Eqns (1) and (2) to eliminate R and 
explicitly introduce the timescale ratio r^. Equating Eqn (7) 
and Eqn (16) then gives 



,(2-/3) 



(17) 

Using /3 ~ 5 from IFerrarese fc MerrittI (|2000l ') then im- 
plies that 



or that 



logM = -log(r^) 



fc, 



(18) 



(19) 



where fc absorbs the constants in the proportionality re- 
lations used to derive Eqn (18). We find that a plot of log M 
against logrr should show a gradient of -1, if the systems 
being considered are indeed quenched by AGN which have 
similar Eddington ratios, regardless of their size. Eqn (18) 
explicitly predicts that, for galaxies which are quenched by 
AGN alone, the quenching efficiency (which is inversely cor- 
related to Vt) increases as the mass of the galaxy increases. 

We now refer back to Figure [S] (top panel), where we 
have plotted the dependence of on galaxy mass for both 
massive and low-mass galaxies. It is clear that low-mass 
galaxies are inconsistent with the energetics of AGN-driven 
quenching, since for these systems Vt correlates positively 
with galaxy mass. Massive galaxies, on the other hand, do 
show an inverse correlation with r^. However, the observed 
gradient in the correlation in Figure [5] is —0.36 ± 0.11, shal- 
lower than the predicted gradient of -1. 

This discrepancy arises due to our assumption, in Eqn 
(13), that the Eddington ratio during the burst phase (n) 
does not change as a function of host galaxy mass. Clearly, if 
jj, decreases, then the rate of energy output (L) also decreases 
and a longer timescale (r2) is required to inject the same 
amount of energy. This, in turn, implies a larger value of r^. 

We now refine our previous analysis by allowing the 
Eddington ratio to vary as a function of mass so that, 

fi^M^ (20) 
With this dependence on M, Eqn (16) becomes 



I1J-/3/3 0/3 



(21) 



„AGN 13 

E oia^.T2- (15) 

Since, the virial theorerrF^I implies that ~ M/ R we 
rewrite this as. 



E (X 



■T2, 



(16) 



The virial theorem states that, for a system with kinetic energy 
K and potential energy U , 2K + U = 0. For a mass distribution 



of mass M and radius R, K • 
that 0-2 ~ MjR. 



and U ■ 



GM 



which implies 



SO that. 



f(37-t-/3~2) ^i^-m 



Using /3 ~ 5 as before then gives. 



M (xr ' 



-1/(1 + 7) 



(22) 



(23) 



Comparison with the observed gradient of —0.36 ± 0.11 
implies that 1 < 7 < 3. We find, therefore, that the ob- 
served gradient in the log A/ vs. logrr correlation can be 
reproduced if the Eddington ratio of the E-l-A galaxies is 
assumed to vary with M. Since the E-fA masses are spread 
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over an order of magnitude, this implies that the Edding- 
ton ratios of the smallest AGN (with host galaxy masses 
M ~ 10^^ Mq) in this sample of galaxies are at least one- 
tenth of those for their most massive counterparts in this 
sample (if 7 ~ 1). 

Finally, it is instructive to calculate if a variable Edding- 
ton ratio (with 1 < 7 < 3) might indeed be expected from 
simple arguments. The luminosity (L) of an AGN is pro- 
portional to the mass accretion rate (M). Spherical, Bondi- 
Hoyle accretion l|Bondi fc Hovlell944 ) implies that M is pro- 
portional to the square of the mass of the central black-hole 
(BH) i.e. 



L (X M (X Mbh- 



(24) 



Since the Eddington luminosity (L^dd) is proportional 
to Mbh we have, 
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The BH mass is typically a consta nt fraction of the mass 
of the bulge (e.g. iHaring fc Rixlliooi '). with 



Me 



Mbulge 

which implies that, 

L 



10" 



(X M, 



(26) 



(27) 



where we have assumed that the E-|-A galaxies in our 
sample are bulge-dominated. Such a scenario implies that 
the expected value of 7 in the analysis presented above is 
~ 1, which is consistent with the derived range (1 < 7 < 3). 

Assuming a more general dependence of M on Mbh 
(M oc Mbh)j would imply that. 



oc M 
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Figure 8. Migration tracks of E+A galaxies in the (NUV — r) 
vs {g — r) colour space (top panel) and the (NUV — r) vs M{z) 
colour-magnitude space (bottom panel). Ages (in Gyrs) along the 
track are shown colour-coded. 



(28) 



A — 1 is equal to 7 in Eqn (20) and the Eddington ratio 
increases with host galaxy mass (as is required to satisfy the 
shallower gradient in the logM vs. logr^ relation) if A > 2. 
For Bondi-Hoyle accretion A ~ 2. 



7 MIGRATION TIME TO THE RED 
SEQUENCE 

Galaxy colours show a pronounced bimodality over a large 
range in redshift. This dichotomy in the colours, com- 
bined with the build-up of the red sequence over time (e.g. 
iBell et al]|2004h . indicates that a central feature of the evo- 
lution of the galaxy population is the net migration of galax- 
ies from the blue cloud onto the red sequence. Therefore, in 
the context of understanding the macroscopic evolution of 
galaxies, it is useful to have an estimate of the typical time 
it takes to complete this migration. 

Given the very recent truncation of star formation in 
these systems, E-l-A galaxies are just about to begin this mi- 
gration process. Assuming there is negligible star formation 



during this phase, we can use our derived SFHs to estimate 
the time it takes for galaxies to move from the blue cloud 
to the red sequence in both the UV and optical colours. The 
use of E-|~A galaxies and the fact that we study the migra- 
tion in both the UV and optical spectrum makes the derived 
migration times quite robust. Note that we define the red 
sequence as the part of the {NUV — r) vs {g — r) colour 
space that hosts the bulk of the red early-type population. 
Comparison to Figure [T] indicates that the red sequence lies 
redward of {NUV - r) ~ 5.5 and {g ~ r) ^ 0.7. 

We estimate the migration times by 'ageing' the best-fit 
model SFHs of each E-l-A galaxy, keeping secondary param- 
eters such as metallicity and E{B — V) constant. We find 
that the migration times vary between 1 and 5 Gyrs (Figure 
[T]). Most galaxies complete their migration within 2 Gyrs 
and the median migration time of galaxies in this sample 
(shown using the solid red line) is ~ 1.5 Gyrs. Finally, in 
Figure |8] we show the migration tracks of the E-l-A sample 
in the {NUV — r) vs {g — r) colour space (top panel) and 
the {NUV ~ r) vs M{z) colour-magnitude space (bottom 
panel). Ages along the track are shown colour-coded. 
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8 CONCLUSIONS 

We have presented the first study of nearby E+A galaxies 
which incorporates their UV photometry. By exploiting the 
sensitivity of the UV to young stars, we have accurately re- 
constructed the recent star formation histories of 38 E+A 
galaxies in the nearby Universe {0 < z < 0.2), by comb- 
ing optical {u,g,r,i,z) and UV data from the SDSS and 
GALEX surveys respectively. 

Wc find that the burst of star formation that domi- 
nates the post-starburst signatures in these galaxies typi- 
cally takes place within a Gyr, which is consistent with the 
presence of A-type stars in these systems. The stellar mass 
fractions formed in this burst are typically high, ranging 
from 20 percent to 60 percent of the mass of the E+A rem- 
nant. The timescale over which this star formation takes 
place is short, ranging between 0.01 and 0.2 Gyrs. The com- 
bination of short timescales and high mass fractions imply 
high SFRs during the burst. We find a tight, positive corre- 
lation between the mass of the E+A galaxy and the implied 
SFR. While low-luminosity {M{z) > -20) E+As have im- 
plied SFRs of less than 50 MQ-yr^^, E+A systems at the 
high-luminosity end {M{z) < —22) exhibit SFRs of greater 
than 300, and as high as 2000 MQyr~'^. The SFRs are com- 
parable to those found in LIRGs and ULIRGs at low rcdshift 
and our results indicate that massive LIRGs at low redshift 
could be the progenitors of massive E+A galaxies like those 
found in our sample. 

We have performed a comprehensive study of the char- 
acteristics of the quenching that truncates the starburst 
in E+A galaxies. In particular, we have studied how the 
quenching efficiency varies as a function of galaxy mass and 
compared the results to scenarios where the quenching is 
due to SN and AGN, which are the typical sources of kinetic 
and thermal feedback invoked in galaxy formation models. 
We have found that in E+A galaxies with masses lesser 
than W^°Mq, quenching becomes less efficient as the galaxy 
mass increases. However, in galaxies with masses greater 
than W^^Mq, this trend is reversed and the quenching effi- 
ciency scales positively with galaxy mass. In terms of M(z) 
and stellar velocity dispersion (u), the reversal occurs at 
M{z) - -21 and a ~ 100 kms"\ 

Since the mass threshold (lO^^M©) where this reversal 
occurs is in excellent agreement with the mass above which 
AGN become significantly more abundant in nearby galax- 
ies, these results can be qualitatively explained in terms 
of AGN and SN being the principal sources of feedback 
that quenches star formation in E+A galaxies above and 
below M ~ 1O^°M0 respectively. In the absence of AGN 
(M < 10^° Mo), the primary source of negative feedback 
arc SN. As galaxies become more massive, the increasing 
depth of the potential well makes it more difficult to eject 
gas from the system, reducing the quenching efficiency. As 
a result, the quenching efficiency shows a negative correla- 
tion with galaxy mass. We have shown that simple ener- 
getic arguments, based on the assumption that quenching 
occurs through mechanical ejection of gas, are able to sat- 
isfy the observed properties of E+A galaxies with masses 
below 10^" M©. This indicates that, in these galaxies, the 
quenching is likely to be driven purely through mechanical 
feedback from SN. 

Once they begin to appear (M > 1O^°M0), AGN be- 



come the dominant source of energetic feedback. Given that 
the AGN luminosity scales with the mass of the black hole 
(Mbh) which, in turn, scales strongly with the central veloc- 
ity dispersion a as Mbh ~ we expect AGN feedback 
to become more cff'cctivc (and hence the quenching efficiency 
to increase) as the galaxy mass increases. We have shown 
that simple energetic arguments, based on the mechanical 
ejection of gas, indeed expect a positive correlation between 
the quenching efficiency and galaxy mass. However, the de- 
rived properties suggest that the quenching efficiencies rise 
an order of magnitude faster with galaxy mass than pre- 
dicted by simple energetic arguments that assume that the 
Eddington ratio does not vary as a function of host galaxy 
mass. However, the correlation between quenching efficiency 
and galaxy mass can be reproduced if the Eddington ratio 
(/j) of the E+A galax;ies is assumed to vary with M as Af' , 
where 1 < 7 < 3. Since the E+A masses are spread over an 
order of magnitude, this implies that the Eddington ratios of 
the smallest AGN (with host galaxy masses M ~ lO^'^M©) 
are roughly one-tenth of those for their most massive coun- 
terparts in this sample (if 7 ~ 1). 

Finally, we have used our E+A sample to estimate the 
time it takes for galaxies to migrate from the blue cloud to 
red sequence. The persistent bimodality in galaxy colours 
over a large range in redshift coupled with the build-up of 
the red sequence over time suggests that the net migration 
of galaxies from the blue cloud to the red sequence is an im- 
portant feature of the macroscopic evolution of the galaxy 
population over time. The use of E+A galaxies, which have 
just truncated their star formation and are on the verge 
of this migrationary transition, together with the use of 
both UV and optical photometry, produces reasonably ro- 
bust estimates of the migration times. The migration times 
are estimated by 'ageing' the best-fit SEDs of each E+A 
galaxy, keeping secondary parameters such as metallicity 
and E(B — V) constant throughout the migration. Wc cal- 
culate migration times between 1 and 5 Gyrs, with a typical 
migration time of ~ 1.5 Gyrs. 

The study of E+A galaxies arguably holds the key to 
understanding many of the processes that shape the evo- 
lution of galaxies. While our analysis is phenomenological 
in nature, this study has attempted to derive quantitative 
insights into some of these processes, especially the charac- 
teristics of feedback that is responsible for modulating and 
truncating star formation in galaxies. While many of the re- 
sults derived here (e.g. timescales over which AGN or SN 
feedback quenches star formation, the comparative efficien- 
cies of the two feedback modes in potential wells of varying 
sizes or the migration times from the blue cloud to the red 
sequence) could prove useful constraints in galaxy formation 
models, future studies of E+A galaxies at high redshift are 
keenly anticipated because they will provide key insights into 
how such processes, that dictate the evolution of the galaxy 
population, evolve over time. 
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